November 25, 2010 



The Hydrodynamic Solution for Flow Profiles in 
a Binary Strong Electrolyte Solution Under an 
External Electric Field 

Byung Chan Eu 
Department of Chemistry, McGill University 
801 Sherbrooke St. West, Montreal, QC H3A 2K6 

Canada 

Kyunil Rah 
IT and Electronic Materials R&D 
LG Chem Research Park 
104-1 Moonji-dong, Yuseong-gu 
Daejeon 305-380, Korea 



Abstract 

In this paper, we follow the general idea of the Onsager- Wilson theory of 
strong binary electrolyte solutions and completely calculate the velocity profile 
of ionic flow by first formally solving the hydrodynamic (Stokes) equation for 
the ionic solutions subjected to an external electric field by a Fourier transform 
method and then explicitly evaluating the formal Fourier transform solutions as 
functions of spatial positions and field strength. Thus the axial and transversal 
components of the velocity and the accompanying noncquilibrium pressure are 
explicitly obtained. They are rare examples for solutions of a hydrodynamic 
equation for flow in an external electric field. The present results make it pos- 
sible to investigate ways to overcome the mathematical difficulty (divergence) 
inherent to the method of evaluating the formal solutions that Wilson used in 
his dissertation on the conductance theory (namely, the Onsager-Wilson the- 
ory) for strong binary electrolytes. Some examples for the velocity profiles are 
numerically computed. They show how ions might move in the ionic solution 
under an external electric field. A possible way to get rid of the divergence caus- 
ing terms for the electrophoretic effect is examined. In the sequels, the results 
obtained in this work will be applied to study ionic conductivity and related 
transport processes in electrolyte solutions. 
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I. Introduction 



The principal initial motivation for study underlying this work was in our 
desire to benefit from Onsager's theories of conductance of electrolyte solutions 
which should be quite relevant to various subjects in the fields of transport 
phenomena in plasmas ) 1 * 2 * 3 ! ionic systems ^ * 5 * 6 * 7 * and semiconductors^^ and 
micro and nano systems ^ 1 * 12 * 13 * consisting of electrolyte solutions and plasmas, 
since the conductivity and mobility of charged species subjected to external 
electromagnetic fields are of great interest regardless of the scales of the system 
size. At quick glance, one might think that Onsager's theories^^on the subject 
matter appear dated at first glance, but the present authors believe there are 
some features that are still as relevant and significant even today as before. In 
pursuing this line of study we have learned that there are things of importance to 
improve on or extend in his theory of electrolyte conductance and new avenues 
to explore as improvement and extension unfold, as will be discussed in the 
course of this line of study. 

In the study of ionic conductivity, one is interested in currents carried by 
ions in solution subjected to an external electric field not necessarily weak. 
As the field strength increases, the conduction current has been observed to 
deviate from the Ohmic law prediction. We in fact observe the well-known Wien 
effect which is a non-Ohmic conductance of electrolyte solutions subjected 
to high external applied electric fiekP^ that defy a linear theory description. 
Currently, such non-Ohmic phenomena are generally observed, for example, in 
semiconductors^ 9 * 10 * 11 * and expected to be operative in micro and nanosystems 
because of the necessarily large field gradients in such systems. Onsager's theory 
is judged to provide valuable insights despite the differences in the sizes and 
states of aggregation of the systems and hence worth a serious study. 

The general ideas underlying conductance in ionic solutions, which originally 
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was Onsager'sf^ consist of two effects: one is the electrophoretic effect and the 
other the relaxation time effect. These can be heuristically and qualitatively 
understood as follows: Let the force on ion j be 

k J=e ,X (j -!,...,«) (1) 

with X denoting the external electric field and ej the charge of ion j. The 
ion of charge ej in the solution will possess an ionic atmosphere of charge — ej, 
and this atmosphere will be subjected to a force of — e^X. This force will tend 
to move the atmosphere in the direction of force — e^X, while the central ion 
j will be carried in the medium in a direction opposite to the motion of ion 
atmosphere. The velocity of this countercurrent may be readily calculated if 
it is assumed that the entire charge —ej of the atmosphere is distributed in a 
spherical shell of radius kT 1 from the central ion and that the motion of this 
sphere is governed by the Stokes k W 19 * 20 ! holding for the motion of a sphere in 
a viscous fluid. Thus, this velocity of the countercurrent is estimated to be 

^ ' k, ej He f 2 ) 

67T?7 67T?7 ' 

where Avj is the velocity of the interior of the shell and 770 is the viscosity of 
the medium. We are thus led to the result that the medium in the interior 
of the shell will be travelling with this velocity, and that the central ion will 
migrate against a current of this magnitude. The deduction of this expression 
qualitatively elucidates the most important part of the effect of electrophoresis. 
One may improve upon it by means of a hydrodynamic method using the Navier- 
Stokes equations. 

The second effect is: if the central ion possessed no atmosphere it would mi- 
grate with a velocity kj/Ci(Gi = friction constant), but owing to its atmosphere, 
the ion is subjected to a force, kj — Akj, where Ak, is the force arising from 
the dissymmetry of the ion atmosphere created by the movement of the ion, and 
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hence it will move with a velocity, relative to its environment, of a magnitude, 
(kj — Akj ) /Q . Consequently, the net velocity Vj of ion j is given by 

(k 7 --Ak,-) 

vj= 1 3 , ,; +Av J . (3) 

For the case of electrical conduction, we obtain 

(e^-X — Ak,-) ej~X.n 



(4) 



Here Akj*^- represents the relaxation time effect. The aforementioned two 
effects underlie a qualitative explanation of the electrophoretic and relaxation 
time effects in conduction. 

The Onsager- Wilson (OW) theory^ is a theory of conductance of strong 
binary electrolyte solutions in an externally applied electric field. It formally 
calculates the ionic conductances by using the Onsager-Fuoss (OF) equations^ 
for distribution functions and the Poisson equation^ 2 ! for the modification of the 
field Ak,- by the presence of the ion atmosphere and, in addition to the afore- 
mentioned two equations, the Navier-Stokes (NS) equation^ — more precisely, 
Stokes equation — for the velocity of the medium to calculate the electrophoretic 
effect. Henceforth the NS equation refers to the Stokes equation in this article. 
For the NS equation the local body force arising from the external electric field 
and many-body effects must be calculated from the statistical mechanics (ki- 
netic theory) of electrolyte solutions. In the OW theory the local body force is 
calculated from the solutions of the OF and Poisson equations. The OF equa- 
tions are coupled differential equations for pair distribution functions of ions in 
the electrolyte solution, which are Fokker-Planck type equations for ions moving 
like Brownian particles in the solution subjected to an electric field. It is useful 
to note that, although the OF equations were originally assumed on physical 
grounds they can be derived by means of the kinetic theory of dense fluids if a 
Brownian motion model is assumed and the many-particle distribution function 
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is suitably approximated.^ 

In his 1936 Yale University dissertation under the supervision of L. Onsager, 
W. S. Wilson^ presented a formal solution for the NS equation obtained by 
means of the method of Fourier transforms with the local body force deter- 
mined from the OF and Poisson equations. The formal solution consists of sine 
and cosine transforms of rather complicated functions. Since the external field 
is assumed directed in the positive x direction of the coordinate system and 
there is an axial symmetry around the x axis, the natural coordinates are cylin- 
drical coordinates (x,p,&), where x is the axis of the cylinder, p is the radial 
coordinate perpendicular to the cylinder axis, and 9 is the azimuthal angle in 
the plane perpendicular to the x axis — parallel to the field direction. See Fig. 
1 for the coordinates defined. The solution therefore is axially symmetric and 
hence does not depend on angle 9. The formal Fourier transform solution for 
the NS equation obtained thereby, however, must be more explicitly calculated 
as a function of x and p to make further progress in theory. Unfortunately, 
the Fourier transform solution is not a kind that can be easily evaluated in 
elementary functions since it consists of rather complicated integrals. Wilson 
considered a special case of position coordinates for evaluation of the Fourier 
transform solution. 

Since the ion j of ion atmosphere is the object of attention and is assumed 
located at the center of the ion atmosphere, which is taken to be at the coor- 
dinate origin, he sets x = and p — in the integrals making up the solution, 
which then become amenable to exact and analytic evaluations. But there is 
one integral among them that clearly gives rise to divergence making the choice 
of x = and p = inappropriate. 

Despite his argument to the effect that it does not contribute to the solution 
because its contour integral vanishes, it is reasonable to say in retrospect that 
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he simply discarded the divergent term and retained only the convergent in- 
tegrals for the calculation of electrophoresis. Despite this troublesome aspect, 
the result of such a calculation was subsequently entrenched in the literature^ 
on conductance of electrolyte solutions, especially, in connection with the Wien 
effect on the ionic conductance. We find that, mathematically, this difficulty 
arises because the some of the integrals involved in the Fourier transform solu- 
tion are not uniformly convergent for all values of the position coordinates. For 
this reason it is not possible to set x = and p = in some of the integrals 
before their integrations are fully performed. 

We, in fact, will find that, when exactly evaluated, some of the integrals give 
rise to functions diverging as x and p approach to the coordinate origin. In other 
words, the solution of the NS equation for the velocity is inherently singular at 
the origin of the coordinate system. It then becomes crucially important to 
re-examine the OW theory to get around this difficulty and suitably revise it 
before making use of it to interpret experiments on conductance, because there 
arises the question of what is really meant by the electrophoretic effect and 
conductance calculated therewith, given a velocity distribution that depends on 
x and p and is manifestly divergent as x, p — > 0. 

In this work, we avoid Wilson's procedure^ of setting x = and p = 
within the integrals before they are evaluated, yet manage to completely eval- 
uate them analytically or reduce them to easily computable quadratures. Thus 
we obtain relatively much simpler velocity formulas than the Fourier transform 
solution that explicitly depend on x and p as well as the field strength, for the 
axial and transversal velocities and the nonequilibrium part of pressure nonlin- 
early depending on the external field strength. As a matter of fact, the results 
obtained for the velocities and pressure appear to be rare examples for the so- 
lutions of the NS equation of a fluid (e.g., electrolytes) in an external electric 
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field when the flow is laminar that is applicable to practical experimental sit- 
uations. The velocity profiles obtained display the manner in which the ions 
move in the applied electric field, and there is much insight to be gained with 
regard to the mode of ionic conduction in electrolyte solutions. In the sequels 
to this work, we will investigate its application to studies of ionic conductance 
and other transport processes in ionic solutions. 

This article is organized as follows. In Sec. II, the formal ( Fourier trans- 
form) solutions of the OF equations for pair distribution functions and the 
Poisson equations for potentials are briefly reviewed to provide the basis for the 
local body force — an electric field in the present case — that is necessary for the 
NS equation for velocity of the electrolyte solution in an external electric field. 
The solutions of the OF and Poisson equations are those of Wilson's in his dis- 
sertation, which have been also presented in a recent tutorial review article^ by 
one of the present authors. In Sec III, the solution of such an NS equation was 
formally solved in the same manner as by Wilson. We repeat his formal solution 
procedure, because his work has never been published, only a brief summary of 
the pertinent results having been given in the monograph of Harned and OwerP^l 
and in the review article by Eckstrom and Schmeltzer.^ His thesis^ itself not 
only skips details of the solution procedures used, making it difficult to follow, 
but also contains some errors, typographical or otherwise, in some important 
steps and results. Moreover, as mentioned earlier, some of the integrals making 
up the formal solution diverge at the origin of the coordinate system, causing 
difficulty to ascertain the true nature and meaning of the electrophoretic effect 
deduced in his dissertation. For this reason, in this work we will solve the NS 
equation afresh, which by itself is quite worthwhile to learn about for its ped- 
agogical value for those not well versed in fluid dynamics of ionic solutions. In 
Sec. IV, the formal solution consisting of a number of Fourier transforms of 
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complicated integrands including Bessel functions^ of complicated arguments 
will be explicitly evaluated by using the method of contour integrations PS We 
thereby obtain explicit velocity profiles as functions of coordinates and field 
strength for axial and transversal flows of the electrolyte solution in a coordi- 
nate system fixed at the center ion of ion atmosphere. Therefore the velocity 
profiles are of the medium with countercharges and relative to the motion of the 
center ion at the coordinate origin; this point should be remembered to avoid 
possible confusion about the motions of ions in the ion atmosphere under the 
external electric field. In the same section, we also calculate the noncquilibrium 
component of pressure which is consistent with the velocity — the solution of the 
NS equation — and depends on the field strength. In this section, we also present 
some numerical examples for the axial velocity profiles that qualitatively display 
how ions move under the external electric field. It reveals an interesting aspect 
of ionic flow in ion atmosphere subjected to an electric field. More detailed 
numerical analysis of the velocity profiles and their applications will be made 
in the sequels where the theory of conductance, transport as well as thermody- 
namic properties, and comparison of the theoretical result with experiment will 
be discussed. In Sec. V, we examine how one might recover the electrophoretic 
effect of the Onsager-Wilson theory from the result of the present work. Sec. 
VI is for discussion and concluding remarks. 

II. Pair Correlation Functions and Potentials 

The fluid system of interest here is a binary strong electrolyte solution sub- 
jected to an external electric field. In the OF theory of transport in electrolyte 
solutions the evolution equations for ionic pair distribution functions are as- 
sumed to obey the OF equations, which are essentially Fokker-Planck type 
equation d 30 * 31 ! in external fields of body force. These equations are coupled to 



8 



the Poisson equations for potentials of interaction between ions. This coupled 
set of second-order differential equations in fact provides an interesting and po- 
tentially productive approach to Coulombic interaction systems subject to an 
external electric field that, in our opinion, has been overlooked in the fields of 
plasmas and semiconductors. 

The OF equations are linear with respect to the potentials because the to- 
tal force involves the assumption that the total field due to the ions and their 
atmospheres may be obtained by linear superposition of separate fields. On- 
sager and WilsorpHMI so i vec j the coupled OF and Poisson equations by means 
of Fourier transforms in a cylindrical coordinate system in which the electric 
field is aligned to the axis of the cylinder. It should be emphasized that lin- 
earizing OF equations with regard to the force does not imply at all that the 
solutions of the coupled set of equations are linear with respect to the field 
strength. The assumption only ensures superposition of fields. In any case, it 
is the point of the Onsager- Wilson theory as far as the OF equations for the 
ion pair distribution functions are concerned. The formal solutions obtained 
are formally exact. They can be employed to obtain the formal solution of the 
hydrodynamic equations for the velocity of the medium. However, the formal 
results are impractical to use for studying transport properties in electrolyte 
solutions, unless they are further evaluated more explicitly. 

Before presenting the solutions it is useful to note the symmetry properties 
of the distribution functions and potentials, which are made use of to construct 
the solutions. The indices j and i will be designated to stand for the positive 
and negative ion, respectively. We denote by r the relative distance vector 
r = Yij = Vi — Yj = — Vji between ions. The ion pair distribution functions are 
then denoted by fji(Tj,Tij) and r^j), where (fc = is the position 

vector of ion k. This implies that for the nonuniform nonequilibrium system of 
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interest the distribution functions are spatially nonuniform with respect to the 
position of ions, e.g., r& in the case of fki(^k, r ik)- The potentials ipj (ry,ry) 
and ipi (ji,rji) obeying the Poisson equations also depend on the position 
of the ion k and the relative distance vector = — ry between ions j and i. 
Henceforth the position vector will be suppressed in the distribution functions 
and potentials for the sake of notational brevity. The Fourier transforms of 
solution ^ 18 ! 21 ! 26 ! for the OF equations and Poisson equations are as follows: the 
distribution functions are given by 

, 2ze 2 n 2 



fii (r) 



irDk B T J 
(l + R) 



da cos (ax) x 



fjj ( r ) = n2 



2R 2 

2ze 2 n 2 
nDk B T J 
(l + R) 



K ° (M + { \J^K (\ 2P ) - l -^-K (X 3 p) 



2R 2 
da cos (ax) x 



2R 2 



Ko (Aip) + (A a p) + l -^-K (X 3 p) 



R 2 



(5) 



(6) 



fij = fji (±r) 

2 2n 2j q'e 



± 



2n 2 r]'ep' 



da cos (ax) 



irDn 2 Jq 
i? 2 



2i? 2 
da sin (ax) ax 

(1-R) 



1 + R) K (X 1P )-^K (X 2P ) 



2R 2 



K (X lP ) + 



R2 K (X 2 p) - -j^Ko (X 3 p) 



(7) 



Here fa (r) and fjj (r) are the pair distribution functions of two identical ions of 
species i and j respectively, and fij (T 1 ") and fji (±r) are the pair distribution 
functions of two different ion species j and i at relative distance — Yj and 



r^, respectively. On the other hand, the potentials of interaction — in fact, 
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the noncquilibrium parts thereof — are given by 
i>j ( ±r ) = -i>i (=F r ) 



2ze r a r \ 
— — / da cos (ax) x 



^ [(1 + R) K (Aip) + (1 - R) K (\ 2P ) - 2 (1 - R 2 ) K (X 3 p)] 

± J da sin (ax) [*sT (Aip) + tfo (A 2 p) - 2Jf (A 3 p)] . 

(8) 

The potentials -0j (±r) and tpi (T r ) are the nonequilibrium potentials beyond 
the equilibrium potentials, which are the Debye-Hiickel potentials.^ 

In the expressions in Eqs. ((5J) — dSJ , x and p are, respectively, the axial and 
radial coordinates of the cylindrical coordinate system; a is the Fourier trans- 
form variable; Kq (Aip) is the Bessel function of second kincP^ with A; defined 
as follows: 

\\ = a 2 + i K 2 (1 + R) , 

\\ = a 2 + 1 -k 2 (1-R), (9) 

A? = a 2 + -k 2 . 
6 2 

Here R denotes 



a = yi-JL_. (io) 

Furthermore, since the external field is assumed to be directed in the positive 
x direction, the functions in the present theory is axially symmetric around the 
x axis, and hence do not depend on the azimuthal angle 9 of the cylindrical 
coordinates. For the binary electrolyte under consideration, if Zj and z% are 
charge numbers, then z = \zj\ = \zi\; e is the unit charge; D is the dielectric 
constant of the medium; p! and rf are defined by 

, zeX , ze 
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with X denoting the held strength, ks the Boltzmann constant, T the absolute 
temperature; k is the Debye parameter of the electrolyte solution 

where n is the density and rife the density of k. For a symmetric binary elec- 
trolyte k is given by 

/ 47T0 6 71 

K = i^T ( Z = N = M) 

The important feature about the ion pair distribution functions and poten- 
tials presented above is their symmetry properties with respect to interchange 
of the ion positions: 

/«(r) = / K (-r), } rj (v) = f n (-v), f ji {±v) = f ij { T v) (13) 

and 

i>j (±r) = ( T r) , (14) 

and the sign ambiguity in Eqs. (O and ((HJ refers to the ion involved, namely, 
positive or negative ion. It should be remembered that the ion positions are 
suppressed in the distribution functions and potentials; for example, /j, (±r) = 
fji{vj,Yij) and ipj (±r) = tpj (rj, ±ry) with r = r,^. It should be reiterated that 
the functions in Eqs. I©-© do not depend on the angle variable 9 owing to 
the cylindrical symmetry of the system under consideration. For the derivations 
of symmetry properties of the ion pair distribution functions and potentials of 
interaction the reader is referred to Refs. 18, 21, and 26. Here we simply note 
that the symmetry properties of Eq. (|13p and Eq. (1141) are for interchange of 
ionic positions or indices for ions. Although they are not necessary for solving 
the hydrodynamic equations, the Fourier transforms given in Eqs. (|5 |) — (|8 |) are 
explicitly evaluated by using the same method for the velocity and pressure in 
Appendix A for completeness. 
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The Fourier transforms on the right hand sides of Eqs. ©-([5]) represent 
nonequilibrium held-dependent parts of distribution functions and potentials. 
Especially, Eq. © must be combined with the Debye potential if one wishes to 
calculate the full potential function of ion j. 

III. Hydrodynamic Equation and Its Solution 

It is well known that ions interacting through long range Coulomb potentials 
produce their ion atmospheres of the mean radius given by the Debye length k _1 
centered around each of them. In an external electric field the ion atmospheres 
of ions interact with the field producing "dressed" local field which in turn 
influences the hydrodynamic flow and movements of ions in the solution. Such 
a dressed body-force is an input when the hydrodynamics of an electrolyte 
solution is sought, subject to suitable boundary conditions. Such a force must 
be calculated to initiate the solution of hydrodynamic equations — namely, the 
NS equation in the present case. It can be calculated by using the potential 
functions given in Eq. ([5]l. In this study we limit the investigation to the 
case of laminar flow, which allows neglecting the nonlinear inertial term in the 
hydrodynamic equation. We assume that there are no body-forces other than 
an electric field applied. 

A. Local Electric Field 

Assume the field is aligned along the x axis. Since the charge density is given 
by the Poisson equation, and the force due to the field on charge density p by 

D X 

F x = pX = -—V 2 i>(r), (15) 
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on substituting the potential function given in the previous section [see Eq. ([5])]. 
we find the local force in the form 

v - />00 -i 

Fx = ~^ J da— {(1 + R)V 2 [cos (ax) K Q (X lP )} 

+ (1-R) V 2 [cos [ax) K (X 2 p)] - 2 (l - R 2 ) V 2 [cos (ax) K (X 3 p)] } 

V" / roc 

-2V 2 [sin(ax)K (X 3 p)}} . (16) 

For a binary electrolyte solution this is in fact the total local force density. 
It is a "dressed" force if we may adopt a modern terminology often used in 
many-body physics. Since the Laplacian operator in the cylindrical coordinates 
(a;, p, 6) chosen is given by 

^dx^ + ],dp' P dp' + ]PW ( ' 

and the gradient operator by 

v -< + *;s + *-5- (18) 

where d x , d p , and Sg are unit vectors in the cylindrical coordinate system, it 
follows that 

/2 / cos(ax) \ , , _ / cos(aa;) \ fl d d , 



V 2 -°r^( K (Xip) = ( Z ^-Pir -o^)K (Xip) ■ (19) 

^ sin(aai) J y lHJ \ say {ax) J \pdp H dp J y H> y ' 

However, because the Bessel function K (Xip) obeys the differential equation^ 

1 d d 



we obtain 
V 



,^- A f)^(V) = 0, (20) 



2 ( Id ) - W - ° 2 ) ( ) «• cm < 21 > 

for Z = 1,2, 3. Therefore the local external force is given by 

F x = Ci cos {ax) K {Xip) + Si sin {ax) K (Xip) , (22) 
14 



where the repeated index implies a sum over the index (Einstein convention) 
and symbols C; and Si are the abbreviations for the following: 

for K (X 1 p) 

for K Q {X 2 p) (23) 
for K {X 3 p) 

and 

for Ko(Xip) 

for K (X 2 p) . (24) 
w (Xl~a 2 ) for K (X 3 p) 

The force F x is an input for the NS equation of interest. The external force zeX 




therefore is dressed up by the extent of the integral factors in Eqs. (|23l) and 
arising from the Brownian motion of ions interacting with the external field and 
among themselves through Coulomb potentials and consequently depending on 
the field and density. 

B. Navier-Stokes Equation 

For a steady flow the Navier-Stokes equation takes the form 

pw ■ Vv - 7?oV 2 v-?7 b V (V • v) = -Vp + F, (25) 

where p is the density, 770 is the shear viscosity, r]b is the bulk viscosity, p is 
the pressure, and F is the body (external) force density. For an incompressible 
fluid V ■ v = and for a fluid undergoing laminar flow of low Reynolds number 
(typically Re = O(10 -6 ) at the field gradient of 1 kVolt/m in aqueous solution) 
the inertial term can be neglected. Thus the Navier-Stokes equations for velocity 
v are given by the pair of equations 

-770 V 2 v = -Vp + F, (26) 
V-v = 0. (27) 

This set is called the Stokes equation by some authors, but we will refer to it 
simply as the NS equation for an incompressible fluid in this work. Note that 
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the presence of an external field makes the pressure nonuniform in space. It is 
interesting to note that if curl of Eq. P2")l is taken, the Vp term vanishes and 
Eq. takes the form 

770V xVxVxv = VxF. (28) 

Since VxV x v = V (V • v ) — V 2 v by vector algebra, the two equations 
and l[2"7|) may be combined into a single equation 

r/ V xVxv = -Vp + F. (29) 

This is equivalent to the Stokes equation, (|2T))) and (|27p for an incompressible 
fluid. For the present problem F —8 x F Xl where 8 X is the unit vector along the 
x axis. 

To solve Eq. (|29p for v, we observe V • v = 0, which means that there exists 
an axial vector A such that v = V x A, where A must depend on position vector 
r and field vector X, both of which are ordinary vectors, that is, polar vectors. 
Therefore A must be a vector function that must also be an axial vector, because 
v is a polar vector. These two conditions are met simultaneously if A = V x a, 
where a = a (r, X) is a polar vector. Note in this regard that the curl of an 
axial vector is a polar vector, and the curl of an axial vector is a polar vector. 
Thus we may write v in a general form 

v = VxVxa + v°, (30) 

where v° is a constant satisfying the appropriate boundary conditions of the 
velocity Note that here a is a polar vector. Since a — > and also v should 
vanish as |r| — > oo, it follows v° = 0. Thus we will set v° = henceforth. 

As the first step to formally solve Eq. substitute Eq. (J3D1) with v° = 

into Eq. (|29p to obtain the equation 

ryoV xVxVxVxa = -Vp + F. (31) 
1G 



To simplify the quadruple curl term on the left hand side we observe that by 
vector algebra 

V x V x a= V (div a) -V 2 a, (32) 
V x V x V (V ■ a) = 0. (33) 

Then take curl x curl of Eq. (|3"2"|) to obtain 

V x V x (V 2 a) = V (V 2 div a) - V 2 (V 2 a) . (34) 

It then follows that 

VxVxVxVxa=VxVxV(V-a) - VxVx (V 2 a) 

= -V x Vx (V 2 a) . (35) 

Upon using Eq. (|3"2"|) in Eq. (|35p and substituting the result into Eq. (|3"Tj) . we 
obtain the equation 

i]o V 2 V 2 a - F = V (r/ V 2 div a-p) , (36) 

which is equivalent to Eq. (|29p or the NS equation. Because the left and right 
of Eq. (f3l)|) are of two different kinds of vectors, the equation is satisfied if 

7 ?0 V 2 V 2 a = F, (37) 
p = po + ?7oV 2 diva. 

In this manner, we have deduced the equation to determine the vector a, namely, 
Eq. (|37p . given the force vector F, Eq. (f22j) in the previous subsection. Thus we 
have obtained the formal solutions for v and p satisfying the NS equation, Eq. 
(|29|) , in terms of vector a. The solution of the NS equation is now reduced to 
that of Eq. (f37|) , a fourth-order differential equation, given F that is provided 
by the OF equations and the Poisson equations. In summary of these results, 
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we have 



v = VxVxa + v= VxVxa, 



(38) 



p = p + 7? V 2 (V • a) . 



(39) 



Vector a is determined by solving Eq. (|3T|) in terms of the local force density 
given by Eq. (|2"2|) . 

Eq. (|38|) differs from Wilson's expression for v in two aspects: a negative 
sign appears in his equation, probably a typo; and the absence of constant 
term vo, which turns out equal to zero. It turns out that the sign error in the 
curl x curl term on the right of Eq. ([38]) was later compensated by another sign 
error in the process of determining vector a described in the following. 

In Eq. (|39[) po is a homogeneous pressure uniform in space, that is, the 
equilibrium pressure. This equilibrium pressure must be either supplied phe- 
nomenologically by using thermodynamics or from the statistical mechanics of 
the electrolyte solution! 32 ! 33 ! For example, it may be calculated from the formula 



where Uij is the intermolecular potential of pair and g®j (r) is the equilib- 
rium pair correlation function. The index i runs over the species in the fluid, 
including ions. Determination of (r) should be made by following the modern 



for p is also absent in Wilson's result for p. 

To solve Eq. (|3"T)l . substitute Eq. (|2"2"j) into the former, which then reads 




(40) 



theory of equilibrium Coulomb (ionic) fluidsP^EU This homogeneous solution 



V 2 (V 2 a) = 8 x Gi cos(ax)K Q (Xip) +S x S t sin (ax)K (\ip), 



(41) 



where 





(42) 
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Recalling Eq. ([2T1) . we find 

cos (ax) K (Xip) sin (ax) K (Xip) 



V z a=d x C r 



S x Si- 



X?-a* 



+ ** ^4*, (43) 



where A* is the homogeneous solution obeying the equation 

V 2 (V 2 A*) = (44) 

with A* = S X A* . The solution A* must fit the boundary conditions at p infinite. 
Thus we choose 



a* cos (ax) K Q (ap) sin (ax) K (ap) 

A = -d x Ci -5 - o x bi- 



Xf-a 2 



(45) 



Therefore we obtain the equation 



V 2 a =S X Q 



cos (ax) [K (Xip) - K Q (ap)} sin (ax) [K (X l p) - K (ap)} 



X 2 ~a 2 



S x Sr 



Xf-a 2 



(46) 



The solution of this inhomogeneous second-order differential equation is sought 
in the form 

cos (ax) sin (ax) 

{h [Ko(Xip) - Ko(ap)] + b 2 [K (j3ip) - K (ap)]} , (47) 



where hi, b 2 , and /3/ are constants determined as follows — note that in Wilson's 
dissertation only the final results for the coeffficients are given. On inserting 
this into Eq. (|4"6"j) we find 

61 (A 2 - a 2 ) K (X lP ) + b 2 (/3 2 - a 2 ) K (/3 lP ) = K a (X lP ) - K (ap) 



[bi (A 2 - a 2 ) - 1] K (X lP ) + b 2 (j3? - a 2 ) K (/3 lP ) + K (ap) = 0. (48) 

Since Kq(Xip), Kq(/3i P ), and Kg(ap) are not equal to zero for all values of p, if 
b\ and b 2 are chosen such that 

1 



(49) 
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and 



lim b 2 (P? - a 2 ) K ((3 lP ) = -K (ap), (50) 

pi— >a 



then Eq. (gJJ is a solution of Eq. (|jg|). Eq. flSOJ) implies 



1 

(ft 2 -a 2 ) 

Finally, we obtain for the solution of Eq. (|4U)) 



h * = -7*T— TV (51) 



a = S x [d cos (ax) + Si sin (ace)] x 



' [K (\ iP )-K (ap)] + M2 * 2^(^)1. (52) 



((A?-a2) 2L UV ir ' " V r " (A 2 -a 2 ) a 2 2 
Here again, there is a sign difference between the last term on the right of 
Eq. (|52|) and the corresponding term in Wilson's equation, Eq. (4.11), in his 
dissertation.^ This particular term gives rise to the aforementioned divergence- 
causing integral, but since it is essentially ignored in his work this sign error 
would not affect his final result for the electrophoretic effect. 

Since the vector a is now obtained, it is possible to calculate the velocity — 
namely, the solution of the NS equation — and pressure by using Eqs. (|38p and 
respectively. For the purpose we calculate diva, V^diva, and V p diva. 
Since for the present problem 

diva= ^ (53) 
ox 

owing to the fact that F =8 X F X and hence a p = ag = identically, we find 
div a = — [Ci sin (ax) — Si cos (ax)] x 

\w^f [KoM " Ko{ap)] + n^] Kl{ap) } ' (54) 

From this follow the expressions 

div a = — [Ci cos (ax) + Si sin (ax)] x 

' [K a (\ip) - K (ap)] + - " P zr Kiifxp) \ , (55) 



(A 2 -«2) 2L uv uv 2(\f-a^Y 
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Vp div a= [Ci sin (ax) —Si cos (era)] ax 
J XiKxiXip) aKi(ap) 



(A? -a 2 ) 2 (A 2 -a 2 ) 2 2 (A; - a 2 ) 



2 /da^A _ 2 r ^ ^ a /aJr (V) 



K (ap) } , (56) 



V &?J = _V [G Sin (a2;) " Sl C ° S \ (A 2 -a 2 ) " 2aR o(ap) j , (57) 
for which we have used the recurrence relation for the Bessel function ;] 28 ! 34 ! 

K' Q (z) = -K 1 (z), 

K[ (z) = -K (z) - ~K X (z) . (58) 

Here the prime denotes the derivative with respect to z. 

C. Formal Solution for the Axial Velocity 

It is now possible to calculate the formal solution for the axial component 
of the velocity. Since 

(V x V x a) x = (diva) - \7 2 a x , 

by using the formulas for V x (div a) and V 2 a x we obtain the formal solution for 
the axial velocity component for all values of x and p: 

v x (x, p, 0) = — Ci cos (ax) x 

A? r , . , ap 



[K (X lP ) - K Q (ap)} + —r—K^ap) 



\(X 2 -a 2 f' UV UV 2(X 2 -a 2 )' 

5; sin (ax) x 

I Xf . , ap 



, 2 [K (X lP ) - K (ap)} + —r—K^ap) . (59) 
(Xf - a 2 ) z ( A i - a ) J 

This formula does not agree with Wilson's Eq. (4.17): Ko(ap) is missing in 

his second term, probably another typo; and the last term is K (ap) instead of 

Ki(ap) as it is here. Since he does not list the sine transform term (the second 

group of the terms in Eq. (|59jl ) it is not possible to compare the present formula 

with his formula for the sine transform part. 
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On substitution of Ci and Si the axial velocity formula can be further sim- 
plified to the form 

zeX , , zeX ( u, \ , , 

v, (x, P , 0) = — g— I c (x, p, + ^ /. x, p, , 60 

where 

/"OO 

I c (x,p, 0) = / da cos (ace) x 
1 



o 



[A?tf (Ai/>) + \ 2 2 K (\ 2P ) - 2X 2 3 (1 - R 2 ) Ko(Mp)} 
^K (ap) + ^K 1 {ap)\ (61) 



and 

I s (x, p, 0) = / <ia sin (ax) 



2a 
^R 2 



Al r^o(Aip) + -^-^ (A 2 p) 



(1 + R) " v ±r/ (1-i?) 
-2AiJf (A3p)] - ^20m K ^ ap) } • (62) 
WilsorPB evaluated the integrals in I c (x, p,0) at x = p — at which point 
I s (x, p, 0) is identically equal to zero owing to the sin (ax) factor. But the last 
integrand in I c (x, p, 0) then gives rise to ^, which yields a divergent result 

f°° 1 

/ —da = oo. (63) 

Wilson argued that this term does not contribute to the electrophoretic effect 
because its contour integral along a semicircle in the complex plane vanishes. 
However, this argument is fallacious because the integrand of this integral does 
not satisfy the Jordan lemma 2 ^'for the contour integral involving a contour along 
an infinite semicircle. Moreover, even if the contour integral along the infinite 
semicircle vanishes, it does not necessarily mean the integral ^da vanishes. 
This, on the contrary, is manifestly divergent and does not need a method of 
contour integration for evaluation. 

To avoid this difficulty we will evaluate the integrals in I c and I s for arbi- 
trary values of x and p. We will then apply the results to the calculation of 
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clectrophoretic effect in the sequel; as a matter of fact, we will have to explore 
a way to obtain it in a finite form. 

This evaluation of v x (x, p, 0) together with evaluations of the transversal 
velocity and pressure for arbitrary values of x and p constitutes the principal 
contribution of this work to the hydrodynamics of strong binary electrolyte 
solutions in an external electric field. The results are the complete solution of 
the NS equation for binary electrolyte solutions in an external electric field; they 
are not only new, but also portend the necessity of a revision of Wilson's result 
for the ionic conductivity IISISI] 

D. Formal Solution for the Transversal Velocity 

By using the relation 

(V x V x a) = V p (div a) - V 2 a p = V p (div a) (64) 

in the case of a p — ag — we find the transversal velocity component in the 
form 

zeX u! zeX 

v P {x, p, 0) = — J c (x, p, 0) - —r— J s {x, p, , (65 

ZTr^r/o k 27M770 

where 

/>oo f 2 

J c (x,p,Q) = da cos (ax) | K 2 R 2 q _ ffl) t Al i l ~ R ) K i(^iP) 

+A 2 (1 + R) K^Xzp) - 2A 3 (1 - R 2 ) K^p)] 

--i^wf^Y (66) 

J s (x,p,0) = J da sin (ax) | ^ 2 ^ 2 [Aiifi(Aip) x 

+X 2 K 1 (X 2 p) - 2A 3 (1 - R 2 ) Ki(X a p)] 
-^- Kl (ap) + ^K (ap)Y (67) 

The integrals in these expressions can be evaluated for all values of x and p in 
a similar manner to the axial velocity, as will be shown. 
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E. Formal Solution for Pressure 

Since for the present system the (nonequilibrium) pressure is given by 

2 / da x \ 

it is easy to calculate it from Eq. (|57|) : 

zeX f°° , a sin (ax) „. „ . 

+ (1-R) K (X 2 p) - 2 (1 - R 2 ) K (X 3 p)] 
daa sin (ace) K$(ap) 

~ t44 r da a2c °* {ax) [K (X lP ) + K (X 2 p) - 2K (X 3 p)] . (68) 

The formula presented above represents a nonequilibrium part of pressure that 
is consistent with the velocity components obtained as the solution of the NS 
equation for a fluid in an external electric field. We have already mentioned on 
how the equilibrium (homogeneous) pressure might be calculate d 32 * 33 ^ see Eq. 
(HOD. 

IV. Evaluation of the Formal Solutions for Velocities and 
Pressure 

The formal Fourier transform solutions for the axial and transversal velocity 
components and pressure, though exact, need explicit evaluation as functions 
of x and p before they can be made to readily indicate their profiles in the 
configuration space (x, p, 0) and applied to study transport and nonequilibrium 
properties of electrolyte solutions. Here we will reduce them to either analytic 
forms or quadratures which can be readily evaluated by means of simple nu- 
merical methods or approximation methods. One can, of course, bypass this 
procedure and simply resort to a numerical computation method to evaluate 
the formal Fourier transform solutions, but this would be rather cumbersome 
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and time consuming computationally, especially because most of the integrands 
involved are singular at some points on the real axis and one has to make use of 
the methods for singular integrals^ to evaluate them. The method of evaluation 
for the integrals discussed below makes it unnecessary to use a direct numerical 
evaluation method. 

It was mentioned earlier that Wilson evaluated the integrals by taking the 
position variables (x, p) at the coordinate origin (x = 0, p = 0) and that such a 
choice of position gives rise to a divergent integral, Eq. (1631) . In the following we 
do not assume a particular set of values for x and p, but evaluate them either 
by means of the method of contour integration or by conventional methods of 
evaluation. 

It is convenient for evaluating the integrals to make use of reduced variables. 
We use the following reduced variables for the purpose: 

9 = (^V^^rio/zeXi^) v, (69) 

t = V2a/n, r — pn/V^, x — xk/V2, £ = ///k, (70) 
V2Xi 



Ul = 



UJ 2 



W 3 = 



K 
K 

V2X 3 



= Vl + * 2 + (71) 



1 + t 2 - y/1 - 2?t\ (72) 



= Vl + t 2 . (73) 



We also use the following abbreviations: 



x/i-^ + v/T+W, S7 3 = v / T 3 y 2 , "= v ^ ■ (74) 
A. Axial Velocity 



The reduced axial velocity then takes the form 

2 V2 



%(x, P ,0) = Ik c + 4=K s , (75) 
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where 



/>oo 

K c {x, r, 0) = / dt cos {tx) 
Jo 

/•OO 

+ 2 / dt cos {tx) 
Jo 

/•OC 

tf s (s, r, 0) = / 
Jo 



(i - 2<e 2 t 2 ) 

[ujlK {u>ir) +LjlK (uj 2 r) - i^LjjKo^r)] 
-w 2 X (rt) + y/fi(rt) 



(76) 



, i sin (ta) 
dt- -^r4rx 



(1 - 2£H 2 ) 

0- + R) 0--R) 

l Z* 00 

~2 / ciiisin(ia;) Ko(rt). 
s Jo 



- 2cj 2 i^o(w3»') 



(77) 



It is convenient to decompose these integrals into various components for the 
sake of separate evaluations as below: 



where 



K c (x, r, 0) = K{ + K c 2 - A^K C 3 - 2K% + rK c 5l 
K s (x, r, 0) - K{ + K s 2 - 2K S 3 - ±K S 4 , 



/>OG 

K{ = dt cos {tx) K (u ir ), 

f°° LO 2 

Kl= \ dtcos{tx)- ^^K {u} 2 r), 

Jo 



"I 


(1 


- 2£ 2 t 2 ) 






(1 


- 2^ 2 t 2 ) 




t 2 LU 2 


(1 


- 2£H 2 ) 



/'OO 

K% = \ dt cos {tx)wlK {rt), 
Jo 

/•OO 

^ = / dt cos (fx) tK^rt), 
Jo 



(78) 
(79) 



(80) 
(81) 
(82) 
(83) 
(84) 
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and 







n 



(1 - 2£ 2 i 2 ) 

/•DO 

K\= \ dtt sin (tx)K (rt). (88) 
Jo 

Henceforth, when the reduced set of variables is used, we will omit the 
overbar from x for notational brevity and understand by x the reduced variable 
x as defined earlier; see Eqs. (lM|) -([T2" j) . 

1. Integrals Amenable to Conventional Methods of Evaluation 

We first treat K%, Kg, and iff, which are all amenable to conventional 
methods of evaluation and give rise to elementary functions of x and p. 

(a) Kl 

The integrals listed above can be evaluated by using the integral representa- 
tion^ of the Bessel function K v {rt) of integer order: 

K v {z) = J°° dse- zcoshs cosh (iss) (\axgz\ < |) . (89) 

On substitution of this integral representation into integral iff and interchang- 
ing the order of integration followed by change of variables we obtain 

7T 7T (2a; 2 - r 2 ) 
K%= — i -j-. 90 

2 ( a .2 +r 2)i/ 2 2(r 2 +x 2 f /2 

(b) Kl 

Upon using the integral representation of if i (rt) and the same procedure as 
for integral iff to evaluate if 5, we obtain 

irr 

Kl = r . (91) 

2 (x 2 + r 2 ) 2 
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(c) Ki 

This integral can be evaluated in the same manner as for K§. We obtain 

KI = ^— T . (92) 

2(x 2 + r 2 Y 

Collecting the integrals evaluated up to this point and expressing explicitly in 
terms of x and r, we find 



2 5 



7T X 

2\/2£( a; 2 + r 2)§ 

- - + ^ + - ^ - ^ (93) 

This contribution of (vi;) acmc to v x represents the fully deterministic part of 
the hydrodynamic velocity that is not associated with the Brownian motion of 
particles giving rise to the dissipative part of the local body force. In fact, some 
of these terms are independent of the field, but (va;) acme is divergent at the ori- 
gin. This divergence, which is an exact result for the integrals partly making up 
?! (x, r, 0), puts Wilson's result for the electrophoresis effect in a rather vexing 
situation, which necessitates careful examination of the flow profiles obtained of 
electrolyte solutions before calculating the electrophoretic effect from them. 

2. Integrals Evaluated by Means of Contour Integration Methods 

There are six integrals whose integrands not only are singular, but also 
involve rather complicated functions and arguments for the Bessel functions 
Ko(Xip) (I — 1, 2, 3) in the expression for the formal Fourier transform solution 
for the velocity. They defy conventional evaluation methods, but, fortunately, 
methods of contour integration may be employed for their evaluation. 

To devise methods of contour integration it is necessary to learn about the 
mathematical properties of the integrals involved. We list their relevant prop- 
erties below: 
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(1) The zeros of the arguments of the Bessel function K„(\ip) for p ^ are 
found to be: 



a — ±iy/ ' k 2 + p! 2 for Ai or, when reduced, t — ±iy2 (1 + £ 2 ) for wi, 

(94) 

a = for A2 or t = for W2, (95) 

a = ±i— ■ p for A3 or t — ±i for W3. (96) 
v2 



The argument of K v {\ip) has branch points at a = ±i-y/ k 2 + /i' 2 , whereas the 
argument of K v {\ip) has branch points at a = and — 00 and the argument of 
K u (X 3 p) has branch points at a — ztin/\/2. 

Thus we may insert a branch cut on the imaginary axis of a plane between 
a = i\J k 2 + p' 2 and a = —i\Jn 2 + p' 2 for the integral of K^(Xip), while a 
branch cut may be inserted along the negative real axis for the integral of 
K v (\2p)i and on the imaginary axis between a = in /y/2 and /a/2 for 

the integral of K„(\ 3 p), respectively. See Figs 2-4 below. 

(2) We recall that Bessel function K v (z) of complex variable z is regu- 
lar in z plane cut along negative real axis! 28 * 34 -! That is, the Bessel function 
is a multi- valued function in the cut plane. In the present case, K (Xip) 
changes discontinuously as the branch cut 



-Ia/k 2 + p ,2 1 iyj k 2 + p' 2 



is crossed, 



whereas K u (\2p) changes discontinuously as the negative real axis is crossed, 



is crossed 



and K u (\3p) changes discontinuously as the branch cut 
on the imaginary axis. Note that the Bessel functions K u (\2p) and K v (ap) are 
defined in a plane cut along the negative real axis. 

(3) We also observe that all the integrands of the singular integrals in I c and 
I s in Eqs. (pjTj) and (f6"2")l are even with respect to a. 

(4) Moreover, for < arg a < ir we find 

- (l-S) <argA / < (\-5) >5 > 0; I = 1, 2, 3) . (97) 
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Therefore in the upper half plane of complex a 



lim K u (Xip)= lim LJLe-^^O. (98) 

foo \ a \— >oo y ZAip 

By this, if contours for the singular integrals are taken along an infinite semicircle 
in the upper a plane, Jordan's lemmeP^l for contour integration along a circle of 
infinite radius is assuredly satisfied. 

(5) Lastly, all the integrands in Eqs. (|8"D )) -([8"2 ]l and (JHSJ) — dHZJ) have simple 
poles at 

a=± $ mt=± ik- (99) 

There is also a branch cut between t = — and t — , but this particular 



branch cut associated with yX — 2£ 2 t 2 does not play a role in the contour 
integrals considered in the present work, because the real axis is not crossed by 
the contours in performing integrations. 

All these properties (l)-(5) together suggest it is possible to evaluate the 
integrals by using methods of contour integration 2 ^! along the closed contours 
of infinite semicircle as depicted in Figs. 2-4. This is indeed a fortunate combi- 
nation of properties of the integrands involved. We will indicate which contour 
applies to which integral at appropriate points in the discussion. Since methods 
of integration will be similar for the integrals involved in Kf and Kf (I = 1,2,3) 
we will illustrate them with the examples of integrals in Kf and K( in Appendix 
A. The results for the rest of integrals will be simply presented with appropriate 
comments on the use of contours in Figs. 2-4 and useful points if needed. And 
then by combining the results for the integrals the final formula for the axial 
and transversal velocities and pressure will be presented. Incidentally, the same 
contours can be used for evaluating the Fourier transforms for the pair distri- 
butions and potentials in Eq. ([S])-© as shown for completeness in Appendix 
A. 
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(d) Kf and Kf 

As prototypes of contour integrals appearing in present work, integrals Kf 
and Kf are explicitly evaluated in Appendix A. Integrals Kf and Kf both have 
simple poles at t = ± (\/2£) \ There is a branch cut along the imaginary axis 
between t = —i\]2 (1 + £ 2 ) and +i\/2 (1 + £ 2 ) and also a branch cut on the 
real axis between t = — 1/\/2£ and t = +l/\/2£, but the latter branch cut plays 
no role in integration since the path of integration in all the contours does not 
cross it. For this reason the latter branch cut is not shown in Figs. 2-4. For 
evaluation of both Kf and Kf the contour in Fig. 2 is used. The results are as 
follows: 



Kf = - 



V2Wl + 2£ 2 ) . ( 
— sin 1 



8£ 3 



K (uJr) 



n f VW+e~) e-*y (1 - y 2 + y/i + 2ey 2 ) 
U dy TTW Io{uJ11 



(100) 



and 



Kf = 71 ^ cos (xt) K (uJr) 



8£ 4 



y(l-y 2 + y/1 + 2^ 2 y 2 



■/ (wir). (101) 



(1 + 2£V) (l + y/l + 2e V 2 ) 

The integrals in these formulas arise from integration along the branch cut. See 
Appendix A for the details of calculation. Therefore 



2 1 V2 



8V2Z 3 

-u 



sin 



V2£ 



cos (xt) K n (ujr) 



K (uJr) 



n rsf*M¥) e-*y (l - y 2 + y/l + 2lV) 
' dy v 



1 + 2£ V 



(l + + 



7 (wir). 



(102) 



(e) and #| 

For evaluation of these integrals the contour in Fig. 3 is used since the inte- 
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grand does not have a branch cut on the imaginary axis, but there is contribu- 
tions from the residues. They give rise to the following results: 



K c 2 =- 



K s 2 = 



V2tt(1 + 2£-i / r 
8?^ Sm 



V2C 



8£ 4 



cos 



V2C 



K (uJr). 



Therefore we obtain 



-K% + -±=K. 



V2 



tt (i + 2e 



sm 



V2£ 



— cos 



(103) 
(104) 

K (Z3r). (105) 



(f) K% and KS. 

The integrands of these integrals involve a branch cut along the imaginary 
axis from z = — i to +i. Therefore the appropriate contour to use is C3 depicted 
in Fig. 4. The results of their evaluation are as follows: 



^ c = -5(l + 2£ 2 ) K Q (Ur) sin ( X 



i6e 



e-'yy 2 (1 - y 2 ) 

2 jo dy rrw fo( " 3r) ' 

^ = ^4(l + 2e 2 )^o(^)cos(-^ 
tt /- 1 e-^y (1 - y 2 ) ( _ 



(106) 



(107) 



Therefore wc find 



-2£ 2 if 3 c - 



(l + 2£ 2 ) X (cJr) 



sm 



cos 



V2C 



n f 1 e-y^y (l - y 2 ) (ygfr-l) 

2i dy TTW /0 



(w 3 r). (108) 
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Summary for the Reduced Axial Velocity 

Collecting the results presented earlier, we obtain the reduced axial velocity 



Vx (x,p,0) = - 



2 i 2N3/ 2 

x + r ) 



+ 



TT (2X 2 - 



2 ( a .2 +r 2)l/2 4 ( a ,2 +r 2)3/2 2 (x 2 +r 2 ) 5/2 

n e-^ (l - y 2 + y/1 + 2£ 2 j/ 2 ) 
dy v ' 



1 - 



1 + 2£V 



- 1 e- a »V2ey (i - y 2 

IJ dy 1 + 2CV 



l + y/i + 2ey 2 ' 

/o(E7 3 r). 



(109) 



We note that the terms made up of trigonometric functions in Eqs. (|102[) , (|105[) . 
and (|108[) cancel each other out. This velocity formula (|109[) is the velocity 
profile of the countercurrent of the ion atmosphere in the coordinate system 
fixed at the center ion of the ion atmosphere pulled by the external electric 
field. The first four terms on the right represent a "deterministic" part of the 
velocity (v x ) and the two integrals involving the Bessel functions Io(uJir) 
and io(u;3r) stem from the Brownian motion part of the local force — namely, the 
dressed-up part of the force arising from the interaction of the ion atmosphere 
and the external electric field. 

B. Transversal Velocity 

On reducing the integrals, we obtain the reduced transversal velocity as 

v p (a:, p, 0) = ^= (J x c + J 2 C - 2J 3 C ) - -A- J 4 C 
l . _ . i 

(110) 
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where the component integrals are defined by 

j; = 1 ^^) {1 _ 2e / )[1+R f ^), (in) 

JS = J Q dtco S (m) {i _ 2e ^ {i _ R) K lM , (112) 
J3=J o dtcos(tx) ^ (j 3 gi(^r), (113) 

pOO 

J% = \ dt cos (fie)tJTi(tr), (114) 
Jo 

and 

/•oo . 

Jf=jf dtsin(fig) ^ ^(^r), (115) 
75 = jf dtsin(fg) ^ g^r), (116) 
■/£ = / dtsin(fg) ^ g!(o;3r), (117) 

p OO 

J 4 S = / dtsin(ts)t 2 i<ri(rt), (118) 
Jo 

/•OO 

J| = / dt sin (tx)tK (rt). (119) 
Jo 

These integrals can be evaluated in the same manner as for the axial velocity 
despite the fact that most of them are given in terms of Bessel functions one 
order higher than Kq (z) appearing in the axial velocity integrals, namely, K\{z). 
The integrals Jf and Jf for I = 1,2,3 are evaluated by using the contours in 
Figs. 2, 3, and 4, respectively, and integrals J|, JJ, and J| are evaluated by the 
conventional method using the integral representation of K\{z). The results are 
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as follows: 



* + 1) K ( x 
= iii wr sin 1 



dy 



e" X!/ y 2 wi/i(wir) 



(1 + 2£V) (l + v/TT2£V) 



J 3 
T c 



(2£ 2 + l) K^uJr) sin 
(2£ 2 + l)K 1 (ZJr) sin 



/\ - xv y 2 (i-v 2 ) T 



a; 

(w 3 r), 



2r 2 (u 2 + l) 5 2(x 2 +r 2 ) 5 
Similarly, we obtain 



J i =^4 (2£ 2 + l)^i(^)cos( 



V 2 (n-« 2 ) 



dy 



1 + 2f V 

^ = ^4(2e 2 + l)^i(-)cos(-^-), 
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3ttx (. 



x — x r — r 
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In summary, we obtain the reduced transversal velocity 
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(124) 
(125) 

(126) 
(127) 
(128) 



xr — r 
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7T rv-v-'w e ^ y yuj\ 
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,2 /- 1 , - x « y 2 (1-y 2 ) r ,_ , 
< X ^ TTw l( " 3r) ' 



ii(wir) 
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As is the case for the axial velocity, the trigonometric function terms in Eqs. 
fT2D j) -([I23 | and Eqs. ([IM j) -([T^5 | cancel in the formula for the transversal 
velocity formula. As is evident from this expression, the transversal velocity 
also diverges at the origin similarly to the axial velocity does. 

C. Pressure 

The expression for nonequilibrium pressure Ap = p — po also can be decom- 
posed into various reduced integrals for the purpose of evaluation: 

A£ = \ (P{ + P^) ~ 2£ 2 P 3 * - PI - -L (pe + p*) + V2CP 3 C , (130) 

where Ap is the reduced pressure 

Ap = M~4^~) (131) 



and 



poo .2 

Pi=J o dt cos (tx) ^^ goM, (132) 

f°° t 2 
Pi = J dt cos (tx) K (u 2 r), (133) 

r°° t 2 

P 3 C = / dt cos (tx) j- t^t^Ko (u 3 r), (134) 
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Pi=J q dtsm(tx) v _ ^ JfoM, (135) 

poo t(\ R) 

P° = J^ dt sin (to) -±—> K Q (u> 2 r), (136) 

poo .3 

Pl= dtsin(tx)- - 75 - 5T Jr ( W3 r), (137) 

Jo 

f-OO 

Pi= dtt sin (tx)K (rt). (138) 
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These integrals are evaluated similarly to the velocities presented earlier: 
Pi = 7=7Z K o(ur) sin 



+ 2/0 * r "TTW /o( " ir) ' (139) 

P 2 = -T^ K °^ r ) sin (-T=z) ' ( 14 °) 



— = — KoCuJr) sin ( _ 

AV2e \V2t; 

+ 1 ['dye-By V 2 
+ 2j Q dye l + 2ey 



dye~ x \ , ;, 2 2 Io(u 3 r), (141) 



p i = TPi ^o(wr) cos 



4£ 2 UV ' VV2^ 



- j ^ ^---L___^ Jo(EJir)i ( 14 2) 

P| = ^ o( ^)cos(-^), (143) 
Pi = -^r^o(wr)cos 



* f 1 1...-™ y 3 



+ U ^TTW 70 ^' (144) 

4 2r3(l + u 2 ) 3/2 2(x 2 +r 2 ) 3/2 
Collecting these results, we obtain the nonequilibrium pressure profile: 

Ap = -j- 

2 ( x 2 +r 2)3/2 

/ •^/2TTTe) e-^j/ (l + V2£y + y/l + 2^y 2 ) 

- - / dy i -5-5 ^/o(wir) 

4i y l + 2^y 2 ° V 7 

+ 2 y rfy — — /o(w3r) ' (146) 

This shows that Ap is also singular at the origin of the coordinates. It is 
significant to observe that the nonequilibrium pressure Ap is generally negative, 
that is, there is a tension that becomes negative infinite at the origin. This means 



37 



the nonequilibrium pressure is compressional. It seems to be a remarkable result, 
probably deserving a deeper consideration. We will report on a further study 
of this nonequilibrium pressure separately. 

By this, we have now shown that the formal solution of the NS equation of a 
binary strong electrolyte in an external electric field can be expressed in terms of 
elementary functions and well-behaved quadratures of regular Bessel functions 
of second kind l v (uJ;r) (I — 1,3) for all values of x and r. These results show 
how the axial and transversal (radial) velocity and the nonequilibrium pressure 
are distributed in the (x, r) space. In other words, they are velocity and pressure 
distributions around the ions flowing in the medium subjected to an external 
electric field of an arbitrary strength. They are hitherto unknown results in 
hydrodynamics of electrolyte solutions. By using these profiles of velocities and 
nonequilibrium pressure we will be able to deduce, in a well-defined manner, 
numerous hydrodynamic consequences of ionic motions in the medium in an 
electric field, subjected to irreversible thermodynamic principles. 

D. Transformation of Velocity to Spherical Coordinates 

In this work, the NS equation has been solved in cylindrical coordinates for 
the reason that the external field has an axial symmetry, and hence the velocity 
and pressure formulas are given in cylindrical coordinates. On the other hand, 
the ions are regarded as cither point charges or hard spheres with spherical 
symmetry. To examine the motion of the medium around the charges (e.g., 
hard spheres) it is convenient to have the velocity expressed in coordinates 
adapted to the symmetry of ions, namely, spherical coordinates. Therefore, it 
is necessary to transform the velocity vector in cylindrical coordinates to that 
expressed in spherical coordinates. This aim can be achieved with appropriate 
transformations for the velocity in cylindrical coordinates (x, r, 9) to that in 
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spherical coordinates (R, i9, ip), which are related to each other by the equations 

6 = <p, 

x = Rcosd, (147) 
r = i?sin , i?, 

where R is the radial coordinate, # is the polar angle, and tp is the azimuthal 
angle of the spherical coordinate system. It should be noted that the roles 
of x and z are switched from the conventional usage. Accompanying these 
transformations, the unit vectors are related to each other as follows: 

S x = cosflSji — siniM,?, 

S r = sin dS R + cos i}6#, (148) 
Se = $<p, 

where 8 X , 8 r , and Sg are unit vectors in cylindrical coordinates whereas Sr, S$, 
and 8 V are unit vectors in spherical coordinates corresponding to coordinates 
i?, and ip. Since the velocity is decomposable in the two coordinate systems 
as 

v = v x S x + v p 5 r + (0) Sg in cylindrical coordinates, (149) 
= vrSr + v$8$ + (0) 8 V in spherical coordinates, (150) 

where vr and v# are velocity components in spherical coordinates, using the 
relations in Eq. (|148l) we obtain the relations between velocity components in 
two coordinate systems: 

Vr = cos^Vj, + sintfvp, (151) 
= — sint?v x +cost?v p . (152) 

Thus multiplying vr with f— 2y/2£/n) 

vr = vr, (153) 
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we obtain the radial component of the velocity in spherical coordinates: 
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i rV 2 ^) e -y Rcos #y 1 - y 2 + Jl + 2£V ' 
9* = - / <2y ^ ■ J -y 



(1 + 2CV 



V2(i + v / i + 2f 




/i(wii?sin i?) 



+ r / rfy- 

/o 



1 e -yRcos'd ni 1 



y 2 (y - V2£) (1 - y 2 ) 
I + 2£2 y 2 /i(w 3 iJsin<?). (156) 



These formulas may be used to calculate the force on the center ion of ion 
atmosphere in the external field. It is closely related to the electrophoretic 
effect we will need for study of conductance. 

In Fig. 5 and Fig. 6, the radial velocity profile vr is presented for two cases 
of £ to show the behavior of the countercurrent in spherical coordinates. In this 
contour map of vr, its magnitude decreases as the color changes from red to blue 
according to the scale shown on the right in the figures. It is, in fact, divergent at 
the origin, meaning an infinite countercurrent, and gradually decreases along the 
axial direction (i? = 0), but as the angle i? increases beyond $ « 40 in the case of, 
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for example, £ = 0.1 it becomes negative when R < 1 as is evident from the right- 
hand corner of Fig. 5. This behavior of the countercurrent does not basically 
change as the field strength £ increases as is evident from Fig. 6 for £ = 3.0. 
The negativity of v# arises from the deterministic part (v) acme defined by Eq. 
(|93[) . which can be negative, but the Brownian motion part of the contribution 
[cos $0 (R, £) — sin #0 t (R, £)] is everywhere positive. This feature is shown 
in Fig. 7 and Fig. 8, in which the contour maps show structures, but do not 
show a negative domain. The color code bar on the right of the figures is for the 
magnitude of the function plotted. In Fig. 9, we show the variation of G (R, 0; £) 
with respect to R and £, showing its decreasing tendency with increasing R and 
also with £ at a given value of R. These figures qualitatively show in which 
region of (R, i?) the countercurrent is least in magnitude, thus offering a least 
resistance to the motion of charges pulled by the external electric field. This 
feature therefore offers a useful insight into studying mobility of ions in the 
electrolyte solutions in the external electric field that, hopefully, can be made 
use of in the light of irreversible thermodynamic principles.. 

V. Relation to the Result by Onsager and Wilson 

Having obtained an axial velocity formula as a function of x and r, it is 
appropriate to see how the axial velocity obtained by Onsager and WilsorP 
for a special position of x = r = 0, namely, the coordinate origin at which the 
center ion of the ion atmosphere is placed, may be recovered. In particular, it is 
important to see in what manner Wilson's formula for the electrophoretic effect, 
/ (£), defined by the relation 

, . z&J£ he , ,. , . 

(0,0,0) = -^— /(£), (157) 

should be understood, given the exact velocity profile formula (|109p . Here 
v x (0, 0, 0) is the axial velocity of ion at the origin obtained by evaluating the 

41 



integrals in the axial velocity formula (15!?)) , or its reduced form in Eq. (f75|) with 
Eqs. ([75)) and (J77J, on taking x = and r = in the integrals. 

At x = r = in Eq. (|59")) , the sine transforms all vanish and the cosine 
transforms contribute well-behaved integrals except for the integral K§ [Eq. 
([HP) ]) which gives rise to a divergent integral, Eq. (p5)) . as pointed out in Subsec. 
III. 4. In Wilson's dissertation,^ provided that the integral K§ is ignored, the 
electrophoretic effect / (£) is found given by 
3 



- (1 + 2£ 2 ) tan" 1 (v^) + (1 + 2£ 2 ) tan" 1 ( y=f j j ■ (158) 

This formula for the electrophoretic effect is in the foundation of the theory of 
Wien effect for binary strong electrolytes.^ In view of the divergent integral 
mentioned, we would like to see in which manner this formula should be obtained 
from the general formula for the axial velocity, Eq. (I109[) , that we have derived 
from the solution of the NS equation. 

The difficulty of getting this formula from Eq. (|109p resides in the term 
related to ( V;c ) acmc [Eq. ®] 

1 r 2 2x 2 — r 2 



(x z + r z ) 



(x 2 +r 2 ) 1/2 2(x 2 +r 2 y (x 2 + r 2 f' 2 

(159) 



which we have termed the "deterministic" part of the axial velocity (— {v x ) acme ) — 
a terminology given to it because it originates from the non-dissipative part of 
the body force in contrast to the other terms in the axial velocity, that is, the 
quadratures in Eq. (I109p . We have multiplied v / 2C/ 7r to this factor because the 
field strength X must be multiplied to obtain the axial velocity. This function 
$ (x, r; £) clearly diverges at the origin of the coordinate system. In other words, 
the axial velocity profiles never passes through the origin at any value of the 
field strength £ because the coordinate origin is a singular point of v x (x, r, 0); 
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in other words, the countercurrent is infinite at the origin. Furthermore, by 
examining the upper bounds of the quadratures, the quadratures in v x expo- 
nentially increase with respect to £ if x is such that xy < Td\r in the interval 



tcrcurrent very large in magnitude in the interval and thus renders the ionic flow 
unrealizable. However, there is a set of trajectories of (x, r) for which xy > uJ\r 
and thus the quadratures are exponentially decreasing as £ increases; these fig- 
ures are not shown in this paper for lack of space, but the velocity profiles are 
instead shown in spherical coordinates in the previous section. Therefore the op- 
timum value of (x, r) must be that of min [v x (x, r; £)] for a given field strength, 
namely, the trajectory of (x, r) corresponding to the valley of the axial velocity 
surface. This trajectory (x, r) x=Xa . r = Tc 6 min [v x (x, r; £)] can be determined as 
a function of £ by a numerical means only. The values of x c and r c are 0(1) in 
reduced units of distance or k _1 in actual units. We have already shown in the 
previous section that this aspect is better displayed if spherical coordinates are 
employed. 

In any case, Wilson's result for / (£) can be recovered, only (1) if the factor 
<£> (ar, r; £) is ignored or the arguments x and r are chosen so that $ (x, r; £) = 
in the definition of electrophoretic effect and (2) if not only x and r are set 
equal to zero in the integrals in Eq. U09\) for v x , but also the terms therein 
that arise from the sine transforms are omitted, that is, the second term in the 
square brackets in Eq. U 60)) given below. We explain it more explicitly: First, 
if we set x = r = in the integrals of Eq. (|109l) wc obtain 



1, yj2 (1 + £ 2 ) in the large £ (field strength) limit. This would make the coun- 



/•V / 2(l+? 2 ) 



(l-y 2 + y/1 + 2£V 
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V2£y 



1 + 2£V 
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l + y/l + 2£V 




(160) 
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For this we note that 

h (z) U=o = !• 

Furthermore, if the second term in the square brackets in the first and second 
integrals in 3 is neglected, that is, the terms n/2£j// ^1 + \f \ + 2£y 2 ^ and 1, 
the resulting integrals are precisely the integrals Wilson^S evaluated to obtain 
the result for / (£) given in Eq. (|f 581) : 



„ 7T 



2 



VW+C) (l - y 2 + v/IT2£V) 



dy~ , , n^9 9 + 2 ^ / d f 



1 „,2 „,2 



r i - r 



2 J -' l + 2£V ' 7 l + 2£V 

(161) 

The integrals here are elementary. It is easy to show that Eyy indeed gives rise 
to /(0 in Eq. (USHD- 

By this, we see that the connection of / (£) to the axial velocity formula in 
Eq. (jl09[) can be made only when we arbitrarily set $ = and also x = and 
r = in the integrals and when some terms therein originating from the sine 
transform terms in the axial velocity are dropped arbitrarily. 

In any case, in view of the fact that the Onsager- Wilson theory for Wien 
effect predicts ionic conductances too small in comparisorP3 with experiment 
it would be useful to examine whether the complete velocity profile formula 
would improve the theoretical prediction. Application to the question of elec- 
trophoretic factor of the axial velocity formula will be made in the sequeP^l 
where conductivity will be studied in comparison with experiment. 

VI. Discussion and Concluding Remarks 

Since ions interact with each other through long range Coulombic interac- 
tions, and ion atmospheres with ions and the external field, the whole body of 
the solution collectively and cooperatively moves subjected to the external field. 
Consequently, the macroscopic behavior of electrolyte solutions in the presence 
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of an external electric field is not simple. 

Their electrical conductivity in particular has attracted a great deal of 
attention in physical chemistry for reasons related to the basic questions in 
chemistry and thermodynamics including the behaviors of electrolytic solutions 
from the early days of physical chemistry. On the basis of Debye's theory of 
electrolytic solutions,^ Lars Onsagerj 14 * 37 ! in particular, elucidated the phys- 
ical mechanisms — the electrophoretic effect and the relaxation time effect — 
underlying the conduction phenomena, which ultimately require solutions of hy- 
drodynamic equations together with appropriate constitutive equations for irre- 
versible nonequilibrium fluxes in terms of thermodynamic forces. The Onsager- 
Wilson (OW) theory of Wien effect,^ and Onsager-Kim (OK) theory^! in the 
case of asymmetric strong electrolytes, depicts the manner in which nonlinear 
field effects on ionic conductivity and mobility can be studied. However, critical 
studies of these theories have not been made in the literature. Even numerical 
studies of, for example, the OK theory, although difficult, have not been made as 
yet. The Onsager-Liu theory^ and the Onsager-Chen theory^ were the most 
recent follow-ups in the aforementioned line, but their works were concerned 
with higher-order density corrections, not the Wien effect itself. 

The ionic association theory approach by Patterson and his collaborators^ 
using the Bjerrum's theory^ of ionic association and Onsager's theory of weak 
electrolytic conduction 3 ^ was based on the electrophoretic effect obtained by 
Wilson despite the question of divergence mentioned earlier. In view of Wil- 
son's electrophoretic effect factor /(£) needing reassessment and probably a 
revision because of the divergence difficulty mentioned, and the notion of ionic 
association that probably can be better addressed by means of the modern sta- 
tistical mechanics of ionic liquids J32133J ^he approach of Patterson et al. would 
require a fresh re-examination. Therefore, the theory of Wien effect and more 
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generally, nonlinear field effects on ionic conductivity of electrolyte solutions is 
not a closed subject as yet despite its long history in the opinion of the present 
authors. 

At this point, it is useful to reconsider, by using the radial velocity compo- 
nent Vfl given in Eq. (|1 54[) , the connection of the present work to the result for 
the electrophoretic effect examined with in cylindrical coordinate in Sec. V. 
Here we would like to consider the radial velocity v# in the direction parallel 
to the external field, that is, d = 0, which then reads 

v fi (R, 0,£) = ^(l-4s)+ 2 ^2 (#> °; • (162) 



R V R 2 . 

Before proceeding further, it is useful to note that if we define a generalized elec- 
trophoretic factor f (R, £) on the basis of the velocity in spherical coordinates 
v (R, ip) = v x cos + v p sin d in analogy to Eq. (|157[) 

then Vfl (i?, 0, £) is related to f (i?, 0, £) by the relation 

v fl (fl,0,0/£ = ^f(fl,0,0. (164) 

If R = y/2, = then 

v fl (V2, 0, e) /€ = ^pf (A 0; e) - 2^29 (y/2, 0; ^) . (165) 

Note that in the actual distance scale R = \[2 corresponds to the radial po- 
sition at 2k _1 , i.e., twice the Debye length, and at this point the factor $ 
[Eq. (I159p ] related to (v x ) , the deterministic part of the velocity, vanishes: 
R=V2,i?=o — 0- And the choice of this particular value of R = y/2 
replaces the step by which the divergent integral was arbitrarily discarded in 
Wilson's procedure. Then the resulting f (V2, 0, £) is rather a similar to /(£), 
which can be shown given by Eyy (Eq. (jTHTJ)). We emphasize that here the 
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choice of R = \/2 is guided by &(R, •&) lij =v /2#=o = 0. We will explore its 
possible irreversible thermodynamics basis in the sequel. 

On the other hand, If 1? = and R = l/\/2, i.e., atv2 times the Debye 
length kT 1 of the radial position, 

v R (1/V2, 0, e) /e = (1/ V2, 0, c) = V2 (29 (1/V2, o; e) - 1) ■ aee) 

These two cases are plotted In Fig. 10 and Fig. 11, respectively, with respect 
to f. We notice that v R (l/\/2, 0,£) /£ >> (\/2, 0,£) /£, or, put in another 
way, the countercurrent is much larger in magnitude in the case of R — 1/ v2 
than J? = a/2- 

These two cases of f (i?,0,£) correspond to / (£) in Wilson's work, but not 
equal to / (£). Evidently, they are well-behaved and finite functions of £, whose 
behavior is reminiscent of that of /(^). Interestingly, f (y2, 0, ^) is closer in 
magnitude to / (£) than f (l/v2, 0, £) . This aspect and the related will be given 
a more detailed consideration in the sequeP^ where application of the results of 
the present work to the theory of conductance will be examined in detail with 
the help of the irreversible thermodynamic principles associated with the flow. 

The present work was born out of our desire to fully comprehend the mean- 
ings of the divergence-causing integral(s) in the formal Fourier transform so- 
lution of the NS equation of binary electrolytes in the presence of an external 
electric field and to find thereby a way to circumvent the divergence in question. 
In this article, to achieve this goal we have fully evaluated the aforementioned 
formal solution of the NS equation without assuming a special case of spatial 
positions, in terms of simple elementary functions plus quadratures of well- 
behaved functions. What we have obtained are the exact velocity and pressure 
profiles in space, which we may apply to study irreversible phenomena in the 
binary electrolyte solutions in the electric field, including electrical conduction 
phenomena. Being a full exact solution without an approximation, it promises 
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to provide a more complete picture of conduction phenomena. This part of the 
study will be made in the sequels to this article. 

The solution [Eqs. (fTDg|) . (ff^l) . (|Ti5|) . and (fT5l|) ] of the NS equation ob- 
tained here is a rare example of rather simple, but explicit solutions for the kind 
of complex, but practical systems, that are amenable to relatively straightfor- 
ward mathematical analysis and computation by a desktop computer, because 
the formulas involved are either elementary functions or comparatively simple 
quadratures free from a singular behavior, such as of poles and other singular- 
ities. They therefore appear to be potentially very useful for gaining insights 
into and describing how ions in electrolyte solutions move in the presence of 
an external electric field. These insights also should help us develop theories 
of electrical conductivity and related transport phenomena in systemJ2~G^ of 
current interest in science and engineering, such as plasmas, semiconductors, 
etc. in electromagnetic fields. 

Acknowledgment 

The present work has been supported in part by the Discovery grants from 
the Natural Sciences and Engineering Research Council of Canada. 



48 



Appendix A. Examples for Contour Integration of Inte- 
grals 

In this Appendix, we use Integrals K\ and K{ [Eqs. ([To]) and ([77]) ] as 
typical examples illustrating the methods of contour integration used in this 
work. Other integrals can be evaluated similarly by using appropriate contours 
given in Fig. 3 and Fig. 4. 

1. Integrals K{ and K{ 

Fig. 2 is used for evaluating this integral. Consider the contour integral 
along contour C\ in complex plane z depicted in Fig. 2 

,2 



( Wl r), (Al) 



where 



1/2 



1 + z 2 + y/l - 2£ 2 Z 2 . (A2) 



Since there is no singularity enclosed by the contour, this contour integral C±Kf 
is clearly equal to zero. Integral CiKf can be decomposed into integrals along 
the paths C_, C+, C, Coo, and along the real axis t. We thus may write it as 



oo 



1 + t 2 + y/l ~ 2^H 2 



1 + z 2 + v/1 - 2| 



2 Z 2 



dz~ 1 _ 2 ^ 2z2 '-e^'Kofa (z) r) 

l + z 2 + v/1 - 2£ 2 z 2 
1 - 2£ 2 z 

(l + Z 2 + y/l - 2£ 2 Z 2 ) 

1 - 2i 2 z 



dz- —5-^5 -e lxz K (uji (z)r) 

dz ~ i t^t^, / -e lxz K (oJi (z)r) 



l + Z 2 + y/l - 2£ 2 Z 2 j 

1 - 2£ 2 z 2 

0. (A3) 



dz- _ 2 ^ 2 '-e^Koiui (z) r) 
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The first integral on the right can be shown to be equal to 2K{ . By the theorem 
of residues^ the integrals C-Kf and C+K^ gives 7rz times residue of (C±Kf). 
Thus we obtain 

C.Kl + C + K{ = ^ ^± 2 ^ sin (-iL) K (Dr) , (A4) 

where 

w = (AS) 

To transform the contour integral C-ftTf around the branch cut along the imagi- 
nary axis we observe that if the phase of the argument of Kq {ui\r) on the right 
hand lip of the cut is chosen equal to zero the phase of the argument on the 
left hand lip is 7ri, so that the argument has the form e m oj\r for the Bessel 
function on the left side of contour C. We now observe that the Bessel function 
can be continued from the left side of the cut to the right side by the following 
continuation formulcP^ 

K {e iw z) = K (z) - nilo (z) , (A6) 

where Iq (z) is the regular solution for the second kind of the Bessel function of 
order 0, Kq(z) being irregular in contrast to Iq (z) being regular. The irregular 
Bessel function K u {z) {v > 0) diverges logarithmically as z —> 0. In series 
representation the Bessel function Iq (z) is given by the formula 

oo / 1 \2m 
m=0 ( m! ) 

This function is finite at z = 0, but it behaves asymptotically as 

I (z)^(2nz)- 1/2 e z [l + O(z- 1 )] (|argz|<). (A8) 
Using formula (|A6|) and changing variable from iy to y, we obtain 



r ^Mi+e) e' x y ( 1 - y 2 + y/l + 2£ 2 y 2 ) 
CK{ = n I dy V 1 + 2gv J -I far) , (A9) 
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where 



1 - y 2 + y/1 + 2^y 2 } V2 • (A10) 



The contour integral along the infinite semicircle Coo vanishes identically be- 
cause Kq {uj\r) vanished 28 * 34 ! as \z\ — > oo. Collecting the results obtained above 
into Eq. (|A3|) . we evaluate the integral Kf in the form: 



V27T (1 + 2e) . ( X , 



^ /-V^ 5 ) e-*" ( 1 - y 2 + v/1 + 2£ 2 y 



d V ' 1 , 0,2 2 ~ J o far) ■ (AH) 



The procedure of evaluating integrals Kf is entirely parallel to the one pre- 
sented above for Kf with the contour in Fig. 2. The result for the reduced 
integral is Kf 

K{ = g£5 (1 + 2£ 2 ) cos (i/V2e) #o (w x (-1/v^f) r)) 

^ /•v / ^ I K 3 ) e-*«» (l - ?/ + + 2£V) 

- ? / dy ^ [ Iq far) . (A12) 

2 Jo (l + 2^V)(l + v / l + 2|V 



The procedures for integrals K%, K£, K^, iff in integrals K c and K s are 
entirely parallel to those presented here for Kf and Kf if the contours in Fig. 
3-4 are made use of. 

2. Evaluation of Pair Distribution functions and Poten- 
tials 

The formal formulas for the pair distribution functions and potentials pre- 
sented in Sec. II can be evaluated by using the same methods as for the velocity 
and pressure. In this part of Appendix A, we present the results for them for 
completeness, although we have not needed them for the purpose of the present 
paper. 
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1. Pair distribution Functions 



With the definition of the symbol 

A/ fe * = (^=y (/« - ™ 2 ) (** = w. y ) ( A13 ) 

we obtain 



A/. = I dy V _____ _ 7q ( - r) 



- f jf ^^ nlV^ o (SJsr) , (A14) 



(1 + 2£V)" 



\/2tt . 
si~ . _ 



A/„ = V sin ( — )^o(S7r) 



+ U dy TTW /o(wir) 



7T r 1 v^ 2 
-2J0 ^(T+W ^ 



and 



y 2 



^^r^irS^) 70 ^- (A16) 



2. Potentials 



With the definition of symbols 



we obtain 



* = * (±r) U^) -^'WEi.) (A17) 



-'I' % Fg?) '' fe,| (' T ^)' (A18) 
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Figure Captions 



Fig. 1 The cylindrical coordinate system employed. 

Fig. 2 Contour C\ for integrals K{ and Kf. This contour also applies to 
integrals Jf and Jf and Pf and Pf . The bold line denotes the branch cut. 

Fig. 3 Contour C 2 for integrals Pf and Pf . This contour also applies to 
integrals Jf and Jf and Pf and Pf . The bold line denotes the branch cut. 

Fig. 4 Contour C3 for integrals Pf and Pf . This contour also applies to 
integrals Jf and Jf and P3 and Pf. The bold line denotes the branch cut, which 
is on the negative real axis. 

Fig. 5 Contour maps of the radial velocity profile for different field strength 
in the case of £ = 0.1. The radial velocity becomes singular at the coordinate 
origin. The magnitude of velocity diminishes from red to blue. 

Fig. 6 Contour maps of the radial velocity profile for different field strength 
in the case of £ = 3.0. The feature of the figure is similar to Fig. 5. 

Fig. 7 The scaled radial velocity profile (P, 1), £) in the case of £ = 0.1. 

Fig. 8 The scaled radial velocity profile (P, £) in the case of £ = 3.0. 

Fig. 9 The cross section of (P, 1?, £) at ■& — plotted against R for 
different values of £. This figure gives an idea of how the (R, £) profiles 
vary with R and the field strength £. 

Fig. 10 The scaled radial velocity f (y/2, 0, £) = at (P = V2, = 0) as a 
function of £. This is equivalent to Wilson's electrophorctic factor / (£), which 
is calculated with v x in cylindrical coordinates instead of spherical coordinates. 

Fig. 11 The scaled radial velocity f (l/\/2, 0, £) = at (P = l/v/2, = 0) as 
a function of £. Comparison of this figure with Fig. 10 gives an idea how the 
electrophoresis might vary with P. 



57 



Figures-II/Figurel.jpg 

Figure 1: 



Figures-II/Figure2 .jpg 

Figure 2: 



58 



Figures-II/Figure3 .jpg 

Figure 3: 



Figures-II/Figure4.jpg 

Figure 4: 



59 



Figures-II/Figure5 . tif 

Figure 5: 



Figures-II/Figure6 . tif 

Figure 6: 



60 



Figures-II/Figurc7.tif 

Figure 7: 



Figures-II/Figure8 . tif 

Figure 8: 



61 



Figures-II/Figure9 . tif 

Figure 9: 



Figures-II/FigurelO.tif 

Figure 10: 



62 



Figures-II/Figurel 1 .tif 

Figure 11: 



63 



This figure "Figurel.jpg" is available in "jpg" format from: 



http://arxiv.org/ps/101 1 .5430vl 



This figure "Figure2.jpg" is available in "jpg" format from: 



http://arxiv.org/ps/101 1 .5430vl 



This figure "Figure3.jpg" is available in "jpg" format from: 



http://arxiv.org/ps/101 1 .5430vl 



This figure "Figure4.jpg" is available in "jpg" format from: 



http://arxiv.org/ps/101 1 .5430vl 



This figure "Figure5.JPG" is available in "JPG" format from: 



http://arxiv.org/ps/101 1 .5430vl 



This figure "Figure6.JPG" is available in "JPG" format from: 



http://arxiv.org/ps/101 1 .5430vl 



This figure "Figure7.JPG" is available in "JPG" format from: 



http://arxiv.org/ps/101 1 .5430vl 



This figure "Figure8.JPG" is available in "JPG" format from: 



http://arxiv.org/ps/101 1 .5430vl 



This figure "Figure9.JPG" is available in "JPG" format from: 



http://arxiv.org/ps/101 1 .5430vl 



This figure "FigurelO.JPG" is available in "JPG" format from: 



http://arxiv.org/ps/101 1 .5430vl 



This figure "Figurell.JPG" is available in "JPG" format from: 



http://arxiv.org/ps/101 1 .5430vl 



